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To elucidate mechanisms of mesoscopic turbulence exhibited by active particles, we experimentally 
study turbulent states of non-living self-propelled particles. We realize an experimental system 
with dense suspensions of asymmetrical colloidal particles (Janus particles) self-propelling on a 
two-dimensional surface under an AC electric field. Velocity fields of the Janus particles in the 
crowded situation can be regarded as a sort of turbulence because it contains many vortices and 
their velocities change abruptly. Correlation functions of their velocity held reveal the coexistence of 
polar alignment and anti-parallel alignment interactions, which is considered to trigger mesoscopic 
turbulence. Probability distributions of local order parameters for polar and nematic orders indicate 
the formation of local clusters with particles moving in the same direction. A broad peak in the 
energy spectrum of the velocity held appears at the spatial scales where the polar alignment and the 
cluster formation are observed. Energy is injected at the particle scale and such conserved quantity 
as energy could be cascading toward the larger clusters. 

PACS numbers: 82.70.Dd, 47.57.jd, 47.63.Gd, 05.65.+b 


I. INTRODUCTION 


Self-propelled particles exist everywhere in nature. 
Their collective behaviors have been studied as subjects 
of nonequilibrium statistical posies Hi- The examples 
include swarming of bacteria ? active turbulence of 
bacteria [Bl, Q , hocking of birds U and schooling of hsh. 
The swimming of such microswimmers as bacteria, al¬ 
gae, and spermatozoa at low Reynolds number has also 
attracted attention from the viewpoint of hydrodynamics 


Some studies have suggested power-law behavior in ac¬ 
tive turbulence, which is mesoscopic turbulence exhibited 
by active particles. Turbulent states of bacterial suspen¬ 
sions have been reported to showpower laws in power 
spectra of their velocity helds p, Q , which is analogous 
to usual huid turbulence. Velocity helds of human crowds 
in a panic are also known to exhibit power laws in such 
quantities as structure function [9| . Power laws in power 
spectra are also reported in numerical studies on tur¬ 
bulent states of self-propelled rods model agent- 

based model pm , and slender rod-like microswimmers 

0 . 

Conditions for the emergence of mesoscopic active tur¬ 
bulence have been investigated by numerical and theo¬ 
retical studies. Wensink et al. 1,0 conducted numer¬ 
ical simulations on self-propelled rods interacting only 
through steric interactions and concluded that a range 
of aspect ratio of self-propelled rods that covers typi¬ 
cal bacterial values can reproduce turbulent states simi¬ 
lar to what was observed in their bacterial experiments, 
which implies that isotropic or spherical particles can¬ 
not exhibit turbulent states. However, in their model. 
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hydrodynamic ehects were dismissed although bacteria 
are swimming very close to each other in huids. GroB- 
mann et al. [ij, [l^ simulated and constructed a kinetic 
theory on an agent-based model with short-range polar 
alignment and anti-parallel alignment at larger distances, 
which also reproduced turbulent states with power laws 
in the power spectra of their velocity fields. Although 
these works suggested sufficient conditions for active tur¬ 
bulence, they omit hydrodynamic viewpoint and thus it 
remains unclear whether such microscopic interactions as 
hydrodynamic interactions in microswimmer systems can 
induce such alignment effects leading to turbulent states. 
It is possible that even spherical microswimmers, whose 
aspect ratio is lower than the values for turbulent states 
reported in the simulations on self-propelled rods model 
[3, [l^ 5 can exhibit turbulent states by realizing appro¬ 
priate alignment effects through any interactions other 
than steric ones. 

Experimental measurements on the power spectra of 
the velocity fields so far have been limited to systems of 
living bacteria, so there is a need to explore power spectra 
in non-living active systems. For the purpose of study¬ 
ing collective behaviors, various experimental systems 
with non-living self-propelled particles have been devised. 
These experiments can help us distinguish between phe¬ 
nomena originating from motility itself and those origi- 
nating from biological activities. Shaken granular rods 
00 or shaken polar disks [HI have been used to 
study collective behaviors, but it is difficult to prepare 
large-scale experiments to remove effects of boundaries. 
Experimental systems with micron-sized particles swim¬ 
ming in fluids have also been invented, although inter¬ 
actions of swimmers are more complicated than those 
in granular systems. Among such microswimmers, self- 
catalytic rods or colloids in hydrogen peroxide solution 
have been the most deeply investigated to demonstrate 
some fascinating characteristics of active Brownian par- 
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tides systems [IM3, although it is difficult to realize 
steady states or vary the motilities of the particles in¬ 
stantaneously in these setups due to the finite chemical 
source of activity. 

Considering these difficulties in other systems, experi¬ 
mental systems which are externally controllable by elec¬ 
tric fields or lasers have been realized and studied |2Q|-[^ . 
One example of such systems is that of Quincke rollers, 
which are symmetrical colloidal particles self-propelling 
under a DC electric field [23-1^. When the particles in a 
hexadecane solution are placed under a DC electric field, 
the symmetry of the charge distribution is spontaneously 
broken, which causes the particles to rotate and move. It 
was found that they show transition from incoherent mo¬ 
tion to macroscopically directed motion. Another exam¬ 
ple is a system of Janus particles, which are asymmetrical 
particles made of dielectric colloids whose hemispheres 
are covered with metal [ 2 ^. When an AC electric field 
is applied to a suspension of the Janus particles, the dif¬ 
ference in the dielectric constant of the two hemispheres 
induces asymmetrical electro-osmotic flow around them, 
which drives the Janus particles to move in a direction 
perpendicular to the electric field [iO, IH, • When the 
suspension is sandwiched between two-dimensional elec¬ 
trodes as depicted in Fig. [U their motions are restricted 
in two-dimensional planes parallel to the electrodes [28| . 

In this paper, we study turbulent states of the Janus 
particles described above. Although experimental sys¬ 
tems of the Janus particles has been realized, experi¬ 
ments with high number density and high voltage, which 
are necessary to observe collective dynamics with high 
motilities, have been challenging because of the techni¬ 
cal difficulty in making the particles to swim smoothly 
without sticking to the electrodes. Hence previous stud¬ 
ies focused on measuring properties of motions of single 
particle or doublets pivoting on an electrode [^, [2814^ . 
Here we improved the surface treatment of the Janus par¬ 
ticles and the electrodes, which enables the particles to 
move at as fast as 61.4 /im/s, which is approximately 
three times faster than the speed achieved by a previous 
study [iO] or typical bacteria such as Escherichia coli and 
Bacillus subtilis. We observed two-dimensional dynam¬ 
ics of the Janus particles at the projected area fraction 
of 0.24. 

We obtained a power spectrum of the velocity field of 
dense suspension of self-propelling Janus particles, which 
has a broad peak. To elucidate the origin of the form of 
the power spectrum, we calculated velocity correlation 
functions and local order parameters for polar and ne¬ 
matic orders. The correlation functions reveal the co¬ 
existence of the polar alignment and the anti-parallel 
alignment interactions, which is analogous to those in 
the agent-based model by GroBmann et al. [nilll. The 
order parameters indicate local cluster formation mov¬ 
ing in the same direction. The range of wave numbers 
where the broad peak is observed corresponds to the spa¬ 
tial scales where the clusters are observed, which means 
that energy is injected and transferred by the particles 


themselves and the clusters in these scales. 


II. EXPERIMENTAL SETUP 


Janus particles used in our experiments are spheri¬ 
cal colloidal particles made of polystyrene whose hemi¬ 
spheres are covered with chromium (Fig. [T]^a)). To fabri¬ 
cate the Janus particles, purchased polystyrene colloidal 
particles (diameter 3 yum. Polyscience, cat^^ 17134) were 
first placed on a glass slide to form a monolayer. Then 
chromium was deposited on this glass slide by thermal 
evaporation to form layers of thickness 20 nm on their 
hemispheres. The Janus particles on the glass slide were 
resuspended in deionized water and washed with soni- 
cation. During the washing procedure, water was once 
replaced with 5 % solution of surfactant Pluronic F-127 
(Invitrogen) to coat the surfaces of the colloidal particles. 
Excess surfactant was removed by replacing the medium 
after sonication and centrifugation. 
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FIG. 1. (Color online) (a) Janus particles used in our ex¬ 
periments are polystyrene beads whose hemispheres are cov¬ 
ered with chromium, (b) Suspension of the Janus particles 
were sandwiched between two horizontal transparent ITO 
electrodes. After sedimentation of the Janus particles, an 
AC electric field was applied in the z direction and the Janus 
particles moved in quasi-two-dimensional horizontal plane (xy 
plane). Images were captured through the bottom electrode 
using a microscope and a CCD camera. 

In our experiments, as shown in Fig. [IJb), suspension 
of Janus particles was sandwiched between two transpar¬ 
ent electrodes (ITO glasses, Sanyo Shinku, Japan). The 
electrodes were separated by stretched Parafilm (thick¬ 
ness: ^ 65 jam). After waiting for a couple of minutes, 
all the Janus particles were sedimented to form a quasi- 
two-dimensional layer near the bottom electrode. When 
an AC electric field is applied, the Janus particles move 
around in this quasi-two-dimensional plane, which is per¬ 
pendicular to the electric field. We observed their mo¬ 
tions through the bottom electrode by using a microscope 
(Olympus, 1X70) with a x40 objective lens (Olympus, 
LUCPLFLN, NA 0.60). Images were captured by a CCD 
camera (Imperx, IPX-VGA210-LMCN) at 200 Hz. 

The origin of the self-propulsive forces of the Janus 
particles is the asymmetrical flow around the particles 
driven by induced charge electro-osmosis (ICEO). Be¬ 
cause no net external forces are exerted on the whole 
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system of the Janus particles and the fluid, total mo¬ 
mentum of the particles and the surrounding fluid have 
to be conserved except in the vicinities of the electrodes. 
Therefore, the counter-action of these asymmetrical flows 
drives the Janus particles to move toward the polystyrene 
hemispheres [^, 0, [^. The motion of the Janus parti¬ 
cles thus caused is called induced charge electrophoresis 
(ICEP). The boundary between the metal side and the 
polystyrene side of each Janus particle is always aligned 
to be parallel to the electric field, because of the hydro- 
dynamic torque exerted by ICEO 0- Hence the self- 
propulsive force is always directed in the two-dimensional 
plane. 

Surface of the bottom electrode was also coated with 
surfactant to avoid sticking. After washing the ITO 
glasses with sonication, they were treated with strong 
plasma cleaning. The cleaned ITO glasses were immersed 
in 5% solution of the Pluronic E-127 for 4 hours and then 
the Pluronic solution was washed away by mildly pouring 
water onto the ITO glasses. 

The gap of the two electrodes was approximately 65 
/im, and the amplitude of applied voltage was 15 V (30 
V peak-to-peak) and its frequency was 1 kHz. Therefore 
the amplitude of the AC electric field was 2 x 10^ V/m, 
which is one order of mag nitude stronger than that in the 
previous studies [^, |^, . This led to higher motili¬ 

ties of the particles, because the self-propelling speed of 
the Janus particles is theoretically and experimentally 
proven to be proportional to the squared amplitude of 
the electric field [20|, [271, l28| . Another reason for choos¬ 
ing this large amplitude was to prevent aggregation in 
this high-density setup, which happened below 10 V at 
this frequency. 

Captured movies were analyzed by detecting all the 
particles. After eliminating noise from the raw images, 
edges of the particles were detected and then binary im¬ 
ages were obtained by filling holes of the edges. In order 
to separate the particles in collision, morphological op¬ 
erations were repeatedly applied to the binary images if 
necessary [32|. The magnification of x40 enabled high 
accuracy of detection, although the observation area was 
smaller compared with lower magnifications. All the par¬ 
ticles were successfully tracked as long as they stay inside 
the frame (Eig. [2]). However, the polarities of the parti¬ 
cles are not detectable with this magnification due to the 
thin layers of chromium. The observation area was 120 
/im X 90/im, which was centered at the sandwiched sus¬ 
pension with the diameter of approximately 5 mm and 
far away enough from the boundaries compared with the 
size of the particles to decrease the effects of boundaries. 
The spatial resolution of images was 0.19 /im/pixel. Du¬ 
ration of the analyzed movie in this paper was 30.7 sec¬ 
onds (6094 frames). Because the self-propulsion speed 
was as high as 61.4 /im/s on average and the particles 
were densely packed, sufficient statistics were obtained 
from these data. 



FIG. 2. (Color online) A trajectory of a Janus particle for 
15 seconds overlaid on the snapshot at the last moment. Due 
to rotational Brownian motion and frequent collisions of the 
particles, their trajectories are not straight but curved. The 
polarities of the particles are invisible with this magnification. 
The scale bar corresponds to 10 fim.. [j^] 


III. RESULTS 

A. Power spectrum of velocity field 

In the previous experimental studies on bacterial tur¬ 
bulence [El, 111 and the numerical studies on a turbulent 
state of self-propelled rods 110, it was suggested that 
power laws in the power spectra might be ubiquitous in 
the velocity fields of the turbulent states of self-propelled 
particles. Such power laws in power spectra are also 
found in fluid turbulence or elastic turbulence [EEj . 
although there are differences in the mechanisms and the 
exponents of these turbulent states. Power-law behavior 
has great significance because it might work as a clue for 
finding scale-free structures or conserved quantities. The 
velocity field of the Janus particles can be regarded as 
a sort of turbulence, because it contains various scales 
of vortices as we see in Appendix. Therefore, it is of 
considerable interest to examine how the power spectrum 
looks like in the Janus particles system. 

The power spectral density E{k) of the two- 
dimensional velocity field of the Janus particles was ob¬ 
tained by applying two-dimensional Eourier transforma¬ 
tion to a velocity correlation function C{R), 

C{R) := {v{t, r) • v{t, r ^ R))t,r , (1) 

/ d^ 7? 

r) ■ v{t, r + R))t,r ,( 2 ) 

where v(t^r) is the velocity of the particle at the posi¬ 
tion r at time t, R is the relative position vector from one 
particle to the other, k is the wave vector in the Eourier 
space, and k is the magnitude of k. All the vectors are 
two-dimensional. The angle brackets mean the average 
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over all the particle pairs at relative position R at every 
time step. The number of calculated pairs was in total 
1.2 X 10^. Because our system is isotropic, the velocity 
correlation function C{R) does not depend on the argu¬ 
ment of R and we can replace C{R) with C{R), where R 
is the magnitude of R (the inset of Fig. [3]). Because the 
distance between the centroids of the particles must be 
larger than their diameters due to exclusive interactions, 
the experimentally obtained velocity correlation function 
has no data points at the distance between 0 jam and 3 
jiim. We interpolated the velocity correlation function 
and obtained the power spectrum according to Eq. © 
(Fig. [3]) im. The wave numbers in Fig. [2] is normalized 
by the wave number 27r/(3/im), which corresponds to 
the particle diameter. The actual calculation procedure 
is summarized in Appendix. [B1 

The obtained power spectrum has a broad peak with 
a lower slope from 0.2 to 0.4 in terms of the nor¬ 
malized wave number, which corresponds to the spa¬ 
tial scales from 7 /im to 15 /im. As the wave number 
goes smaller, the slope of the power spectrum gradually 
changes roughly from the exponent 0.5 to 1. The ex¬ 
ponent 1 in the smaller wave number regime is trivial 
because the Fourier transform of the correlation function 
cannot detect large wave length fluctuations sufficiently 
and the integrals in Eq. m and Eq. (jB7|) behave as con¬ 
stants with respect to /c, which results in E{k) oc k. In 
bacterial turbulence [5|, the power spectrum with the 
exponent 5/3 was suggested, but it is inappropriate to 
compare the forms of their power spectra precisely due to 
the limited observed ranges of the wave numbers for both 
cases. Furthermore, in the study on agent-based model 
with short-range polar alignment and anti-parallel align¬ 
ment at larger distances by GroBmann et al. [a El, it 
was demonstrated that the scalings in the power spectra 
depend on the choice of interaction parameters and thus 
the exponents and the form of the power spectra are not 
universal. 


B. Velocity correlation in detail 

We look into the velocity correlation in detail to inves¬ 
tigate the origins of the form of the power spectrum. In 
the inset of Fig. O the oscillation of the velocity correla¬ 
tion function C{R) is observed, indicating the existence 
of structures with finite length scales. This form of C{R) 
is intriguing because GroBmann et al. also observed quite 
similar oscillation in C{R) of the turbulent states of their 
model ini- To quantify the angular information, which 
was dismissed in C{R), we define another velocity corre¬ 
lation function 

C'{R, Ip) := {v{t, r) ■ v{t, r + -Re^))^,,, (3) 

where is a unit vector in the direction of ^/-harg t;(t, r). 
C'{R^'ip) can extract the information on how the other 
particles are moving from a viewpoint of a reference par¬ 
ticle. 



FIG. 3. (Color online) Log-log plot of the power spectral 
density of the velocity field of the Janus particles. The power 
spectrum was calculated according to Eq. © using the ex¬ 
perimentally obtained velocity correlation function defined as 
Eq. O. The velocity correlation function C{R) is shown in 
the inset. The blue and green lines are slopes with the expo¬ 
nents 0.5 and 1 respectively to guide the eye. Wave number is 
normalized by the wave number corresponding to the diame¬ 
ter of the Janus particle. The vertical axis is not normalized. 


The velocity correlation function C'{R^'ip) shown in 
Fig. m indicates that anti-parallel alignment in the trans¬ 
verse directions and polar alignment in the longitudi¬ 
nal directions coexist. The negative correlation in y' 
( 7 / = ±7r/2) direction, which is perpendicular to the di¬ 
rection of motion, indicates that the particles next to 
each other tend to move in the opposite directions. This 
anti-parallel alignment can be partially explained by con¬ 
sidering hydrodynamic interactions between the Janus 
particles in a sense that the parallel configuration of two 
Janus particles moving in the same direction is unstable 
(see Appendix. Ofoi* a detailed explanation). The pos¬ 
itive correlation in x' direction, which is parallel to the 
direction of motion, indicates that the particles are fol¬ 
lowing other particles in front of them. In the study by 
GroBmann et al. 0E1 , it was shown that the coexis¬ 
tence of polar and anti-parallel alignment lead to effec¬ 
tive negative viscosity and can trigger active turbulence. 
This coexistence might cause effective negative viscosity 
in suspensions of the Janus particles and account for the 
turbulent state in the Janus particles system, although 
the angular dependence of the alignment interactions is 
different from that of their agent-based model. 


C. Local order parameters 

In order to quantify structures and orders emerging 
from the two-body velocity correlations which we have 
seen in the previous subsection, we calculate two locally 
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(a)Colormap of the velocity correlation function 
C'iR.'ip) 



(b)Velocity correlation in x' or y' direction 


FIG. 4. (Color online) (a) Colormap of the velocity correla¬ 
tion function C'{R,'ip) defined in Eq. (|3]). The coordinate is 
x' = Rcos^jj and y' = i?sin^. Note that x' and y' are differ¬ 
ent from the axes defined in Fig.[T] The filled gray circle rep¬ 
resents a reference Janus particle, whose direction of motion 
(the red arrow inside) is rearranged to the right (x' direction). 
The dashed circle corresponds to twice the radius, which the 
other particle cannot enter due to exclusive interaction. In re¬ 
ality, repulsive electrostatic interactions, surfactants absorbed 
on the colloidal surface, and the suspending fluid dominated 
by Stokesian dynamics keep the particles apart and they can¬ 
not get close to this dashed circle. No particles are detected 
in the white region, (b) Slices of (a) at y' = 0 {x' direction, 
parallel to the velocity, red circles) and x' = 0 (y' direction, 
perpendicular to the velocity, blue squares). Negative correla¬ 
tion in y' direction shows anti-parallel alignment and positive 
correlation in x' direction indicates that the particles are fol¬ 
lowing other particles in front of them. 


defined order parameters: an orientational order parame¬ 
ter I (e^^) I and a nematic order parameter | (e^^^) |, where 0 
represents the direction of instantaneous velocity of each 
Janus particle (Fig.[5|). These order parameters can clar¬ 
ify whether nematic ordering or lane formation is ob¬ 
served as they are reported in numerical studies with 
anti-parallel alignment interactions HISISl. We calcu¬ 
late these order parameters locally in a L x L grid to make 


histograms of their values. The grid size L is varied from 
9 /im to 15 jivd. The obtained probability distributions 
are shown in Fig. [6l /(r) plotted in Fig. [6] are probabil¬ 
ity distribution functions (pdf) with respect to moduli of 
(e^^) and (e^^^), 


1 

/(|(e'^‘^)|=r)oc—d,/.pdf(|(e-^) 


= r, 


(4) 


arg(e-^) =(/.), 


where /i: = 1 or 2 for the orientational order parameter 
or the nematic order parameter respectively, and arg is 
a function returning an argument of a complex number. 

As the grid size is decreased, non-zero peaks gradu¬ 
ally appear and grow higher for both order parameters. 
This indicates that the particles are locally moving in the 
same direction forming structures like trains but the op¬ 
posing lanes are not formed. The particles exhibit local 
orientational order, although they do not have global ori¬ 
entational order, nematic order, or positional order (see 
Appendix. p. When we examine the captured movies 
carefully, some local clusters moving in the same direc¬ 
tion for a while can be found frequently ( Fig. [7| and 
supplemental movies [ssl, [39|). The fact that the non¬ 
zero peaks for the two order parameters are observed in 
various grid sizes suggests that there exist various scales 
of these clusters up to 12 /rm in size. These structures 
are formed and broken up at all times, changing their 
cluster sizes. 

The length scales where the broad peak is observed in 
Fig. [3] corresponds to those where the non-zero peaks are 
observed for the two order parameters in Fig. [6l Fur¬ 
thermore, the strongest positive correlation in the longi¬ 
tudinal direction in Fig. SJat x = —7 /im corresponds to 
the highest position of the broad peak in the power spec¬ 
trum. Thus we conclude that the formation of clusters 
of particles moving in the same direction is a possible 
cause of the origin of the broad peak found in the power 
spectrum of the velocity field of the Janus particles. 
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FIG. 5. (Golor online) The local average values of the orienta¬ 
tional order parameter |(e^^)| and the nematic order parame¬ 
ter |(e^^^)| can distinguish local states of collective behaviors. 
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(b)Probability distribution of |{e^^^)| 

FIG. 6. (Color online) The probability distributions of the 
order parameters dehned as Eq. o are calculated at the scale 
of L = 9,10,12,15 fim. respectively. When L is smaller than 
12 /xm, both order parameters have non-zero peaks meaning 
that the particles are locally moving in the same direction but 
opposing lanes are not formed. These peaks gradually grow 
as L is decreased. In the calculation of order parameters, 
the regions in the grid including fewer than three particles 
are neglected in order to extract the results coming from the 
collective dynamics. 


IV. DISCUSSION 

In the usual fluid turbulence, according to Kol¬ 
mogorov’s law the exponent of the power law ob¬ 
served in the power spectrum is negative, while the slope 
in our experiment on the Janus particles is positive. The 
crucial difference is in the spatial scales where energy is 
injected into the systems. In the usual fluid turbulence, 
energy is continuously injected externally by boundary 
conditions such as pressure gradient and shear stress. 
Consequently, the energy input occurs at the scale of its 
system size, which corresponds to the smallest wave num¬ 
ber in the spectrum, and this energy is transferred to the 
larger wave number regime, finally dissipating into heat 
at small eddies. In contrast, in the active turbulence of 



FIG. 7. (Golor online) (a) In a snapshot of the experimental 
movie, each particle is colored depending on its direction of 
instantaneous motion. Spaces between the particles are also 
colored according to average velocities of the nearby particles 
with Gaussian weight (see Appendix. Ofor detail). We can see 
some regions connected with same colors [s^. (b)Trajectories 
of the particles moving in the same direction are overlaid on 
the snapshot at the last moment. Different colors indicate 
different clusters [s^. The scale bars correspond to 10 /xm. 

the Janus particles or bacteria (Ef, energy is injected in 
the bulk by the particles themselves, whose length scale 
is the smallest in the system. This means that energy in¬ 
put is occurring at the largest wave number in the power 
spectrum. According to our experimental results show¬ 
ing that the Janus particles are forming clusters with the 
particles moving in the same directions, it is possible that 
some conserved quantities, which are energy and enstro- 
phy in the case of fluid turbulence, are cascading from 
the particle scale to the larger structures. 

Since the beginning of studies on active matter, most 
effort has been put on describing macroscopically coher¬ 
ent states or transitions between ordered states and dis¬ 
ordered states. However, to the best of our knowledge, 
there are no experimental systems of microswimmers ex¬ 
cept for the Quincke rollers that exhibit macroscopically 
directed ordered phase, which is often obtained in many 
numerical studies including Vicsek model mu. It is 
of great importance to understand from the microscopic 
point of view why it is difficult for microswimmers to 
align in the same direction. 

The interactions dominating the collective dynamics 
of self-propelled colloidal particles under an electric field 
are both hydrodynamic and electrostatic ones. In the 
case of the Quincke rollers, their electrostatic interaction 
had been already understood, so the calculations con¬ 
sidering the hydrodynamic interaction clearly explained 
the mechanism of the transition to the coherent motion 
[ 2 ^. However, as is mostly the case for theoretical cal¬ 
culations, their calculation on the hydrodynamic inter¬ 
actions relies on far-held approximation, in which each 
microswimmer is represented as a point singularity 0. 

Near-held hydrodynamic interactions are much more 
important to understand collective behaviors of swim¬ 
mers than far-held ones, especially in the case of spher- 
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ical particles moving very close to a surface such as our 
Janus particles. Because spherical particles do not have 
strong steric alignment interactions, hydrodynamic in¬ 
teractions dominate their dynamics. When the particles 
are moving very close to a surface or confined between 
two walls, far-held how is strongly suppressed and only 
near-held interactions play an important role in their dy¬ 
namics mLin. When we neglect the near-held inter¬ 
actions and the electrostatic interactions, the remaining 
far-held interactions cannot account for the diherence in 
the macroscopic behaviors of the Janus particles and the 
Quincke rollers. Their far-held interactions are the same 
because the two electrodes close to the particles damp all 
the singularities of how except source dipoles created by 
each particle [10, |42|. Therefore, the diherence in their 
macroscopic states should be explained in terms of their 
near-held hydrodynamic interactions or electrostatic in¬ 
teractions, which must be very diherent between the two 
due to the diherence in the ways how their symmetries 
are broken and the direction of self-propulsion is deter¬ 
mined. 

Because of the difficulties in the theoretical calcula¬ 
tions on near-held hydrodynamics, numerical simulations 
have been performed to understand the hydrodynamics 
of swimmers [i^, In these studies, microswimmers 
are represented as squirmers, which have spherical bod¬ 
ies with slip boundary conditions [i^. This simplihed 
model is useful for studying our Janus particles because 
their slip velocity is theoretically calculated in the low 
frequency limit (see Appendix. El- Numerical simu¬ 
lations on Janus squirmers with electrostatic interactions 
are now in progress [46|. 


V. SUMMARY 

To study turbulent states of active particles by us¬ 
ing non-living microswimmers, we realized an experimen¬ 
tal system of Janus particles propelling in a quasi-two- 
dimensional plane with high density and high motilities. 
The area fraction of the particles can be set as high as 
0.24, which has not been achieved because of aggrega¬ 
tion of the particles in the previous studies [28|. Their 
average self-propelling speed reaches 61.4 /im/s, which is 
about three times faster than usual bacteria or the Janus 
particles in the previous studies [lO, [HI, lUl • We tracked 
all the particles in the observed area and calculated the 
power spectrum of their velocity field which shows the 
broad peak ranging over a length scale of several times 
the particle size. 

Calculation of the velocity correlation functions re¬ 
vealed the coexistence of the polar alignment and the 
anti-parallel alignment, which had been demonstrated to 
be responsible for the effective negative viscosity of ac¬ 
tive suspensions and the emergence of active turbulence 
[nllil. However, it remains to be elucidated how these 
alignment effects arise from the microscopic interactions 
such as the hydrodynamic and the electrostatic interac¬ 


tions. 

The probability distributions of the order parameters 
indicated the formation of local clusters moving in the 
same direction. Judging from the correspondence be¬ 
tween the spatial scales where the broad peak is observed 
and those where the clusters are found, the cluster for¬ 
mation can be understood as one of the possible causes 
of the power law. 

In future, we shall focus on measurements of the prop¬ 
erties of their dynamics by changing control parameters 
such as the amplitude and the frequency of the applied 
AC electric field, and the number density of the particles. 
In this system the particles tend to aggregate at smaller 
amplitude of applied voltage, so we have not successfully 
analyzed how the collective dynamics of our Janus par¬ 
ticles depend on the amplitude and the frequency of the 
electric field. Because we can tune the interactions, the 
self-propulsion speed, and the direction of motion of the 
Janus particles by controlling the amplitude and the fre¬ 
quency of the electric field , it shall serve as an ideal 
experimental setup to figure out some universal laws of 
collective motion of self-propelled particles from their mi¬ 
croscopic interactions. 
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Appendix A: Vorticity map 

In order to examine how turbulent the velocity field 
of the Janus particles is, we make a color map of the 
vorticity of the velocity field (Eig.[8j). The velocity of each 
particle is calculated from its trajectory, and then the 
time series of velocities are smoothed by taking moving 
average over 21 frames. Because the velocity is originally 
defined only on the centroids of the particles, we apply 
Gaussian filter to obtain the velocity field. The standard 
deviation of the applied Gaussian filter is chosen to be the 
radius of the Janus particles. Various scales of vortices 
with pos itive and negative vorticity coexist and evolve in 
time 113 . 


Appendix C: Flow field and hydrodynamic 
interaction 
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FIG. 8. (Color online) Color map of the vorticity. Various 
scales of vortices with positive and negative vorticity coexist. 


Appendix B: Calculation of Spectrum 


When we consider an isotropic two-dimensional sys¬ 
tem, the velocity correlation function defined as Eq. m 
is independent of the angle of R, leading to 


C{R) = {v{t, r) • v{t, r + R))t,r‘ (Bl) 


Calculation of its power spectral density defined in 
Eq. @ is the following: 


Eik) = 


k 

27r 


/ 

L 


CX3 p2tt 

RdR / d(j> 
Jo 


(B2) 

(B3) 


where 0 is the argument of R. We can calculate the part 
of the integrand with 0 by using the periodicity of cosine 
function, 


/ d</> Q-ikR cos 

Jo Jo 

p27T 

= / d</> (B5) 

Jo 


= 27rJo{kR), 


(B6) 


where a is the argument of k and Jo is the 0th order 
Bessel function of the first kind. By substituting Eq. (jB6p 
into Eq. (|B3p . we obtain the final result, 

POO 

E{k) = k dR C{R) • R Jo{kR), (B7) 
Jo 

We can calculate the power spectrum of an isotropic 
velocity field from its correlation function by using 
Eq. 


Janus particles are interacting with each other through 
hydrodynamic and electrostatic interactions. The elec¬ 
trostatic interaction is quite complicated, because elec¬ 
trical double layers around the particles have their own 
characteristic charging time and thus their thickness de¬ 
pends on the frequency of the applied AC electric field. 
On the other hand, the hydrodynamic interaction can be 
understood by solving the Stokes equation. Erom the hy¬ 
drodynamic interactions, we can conclude that the par¬ 
allel configuration of two Janus particles moving in the 
same direction is unstable. 

We can solve the Stokes equation to obtain the flow 
field around a Janus particle by assuming slip velocity 
on the particle, if we neglect the nearby electrodes and 
only consider one Janus particle suspending in fluid. The 
flow around the Janus particle is dominated by Stoke¬ 
sian dynamics, because the Reynolds number of the flow 
is as low as 10“^. On the surface of the Janus parti¬ 
cle, the fluid is driven by ICEO. The ICEO flow on the 
polystyrene side is considerably weaker than that on the 
metal side, so it i^ustifled to neglect the flow on the 
polystyrene side [^. We regard that this ICEO flow 
exists only in the vicinity of the surface, which is the 
thin electric double layer with the thickness of its Debye 
length [^. We replace the non-slip boundary condition 
on the metal hemisphere to the slip boundary condition 
with the slip velocity calculated in a theoretical study 

na. 



FIG. 9. (Color online) The spherical coordinate used for cal¬ 
culation. The origin is at the center of the Janus particle. 
The right hemisphere {x > 0, or |^| > 7r/2) is covered with 
metal. Note that the axes defined in this Appendix. [O are 
different from those used in Figs. [IISIEI 


The method to solve the equation is the following ji^ . 
First of all, we consider the particle frame, in which the 
particle is at rest. We take a spherical coordinate as 
depicted in Fig. [9l The slip velocity on the surface of the 
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Janus particle Ug is given as 

Us{ 0 ,(f>) = 

( —^Uo{sm29cos^ (f> 0 — sin9 sm2(f> (|0| < 7r/2) 

\0 (l^l>^/2), 

(Cl) 


where Uq is the typical speed of ICEO flow, and 6 and 0 
are the unit vectors with +6> and +0 directions respec¬ 
tively. In order to solve the Stokes equation, we need 
another boundary condition, which is the velocity at in¬ 
finity Uoo- Because we are in the co-moving frame with 
the particle, the velocity at infinity has the opposite sign 
and the same amplitude of the velocity of the particle 
t/iCEP in the laboratory frame. C/icep can be calculated 
from the slip velocity and it leads to 

If 9 . 

'^oo = -UicEP ^ ^ = —Uox, (C2) 

where t/jcEP is the velocity of the particle in the labo¬ 
ratory frame, dft is an element of solid angle, and x is 
a unit vector in the direction of x axis. We expand the 
boundary conditions on the surface with spherical har¬ 
monics 0), 



-rW ■ {U - Uoo)\r=a 

r • V X 


oo I 

Y, E ^i,mYr{o,4>) (C3) 

1=0 m=—l 
oo I 

E E (C4) 

1=0 m= — l 
oo I 

E E ci^mYr^oA) (C5) 

1=0 m= — l 


where r is a position vector in the spherical coordinate, 
r is the magnitude of r, a is the radius of the Janus par¬ 
ticle, (j)) is a spherical harmonic function of degree 

I and order m, and and Ci^rn are coefficients 

of expansion. By using Eq. m, Eq. (jC2[) and an or¬ 
thogonality of spherical harmonics, we can obtain the 
coefficients as 


^27r ^TT 

Ai,m = d(p d9 
Jo Jo 

Jo Jo 


d9 


-•(m-Woo) Yri9,cl>) 

T r=a- 

(C6) 

-^f/ocos0 Yr{9,4>) (C7) 


p27r pTT 

Bi,^= d4> d9[-rv ■{u-Uoo)l^jYr{9,4>) 

Jo Jo 

(C8) 

= J J d6> I ^^0 [(2 cos^ ^ — sin^ 6>) cos^ 0 

-cos 2(1)]} Yi^ {0,(1)) 
(C9) 


p27T 

on 


■■ / dct> 

/ d6»[r • V X (M-Moo)lr=J 

YnoA) 

Jo 

lo 

(CIO) 

p2TT 

r "2 


■■ / d(j) 

/ d9-0-Yr{9,4>) 

(Cll) 

Jo 

to 


0. 


(C12) 


Ai^rn has a non-zero value only at Ai^q = — ^ 
Absolute values of non-zero Bi^rn are plotted in Eig. 
which decay algebraically with the increase of 1. 



Degree I 


FIG. 10. (Color online) Values of the coefficients of spherical 
harmonics (a) Bi^rn in a linear plot and (b) \Bi^rn\ in a log-log 
plot. Bi^rn has non-zero values only when m = 0, ±2. For any 
/, = Bi ^-2 holds due to symmetry in the sign of J). 


Because the expansions are infinite series, we take up to 
I = 50 terms to calculate the flow field. Judging from the 
values of ^ 2,0 = 1-78 and ^ 49^0 = —0.0456, truncation 
error can be estimated at approximately less than 3 %. 
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By defining six kinds of harmonics, 


i 


x{i,e,4>) 

:= ^ Ai,^Yr 

m=—l 


(C13) 

Y{i,d,4>) 

1 

:= ^ Bi,mYr 

m=—l 

{OA) 

(C14) 


1 

:= ^ Ci,mYr 

(^,0) 

(C15) 


_ 7?(2n- 1) 1 / 

a\”+i 


71 + 1 a \ 

-r) 



X [(71 + 2)A(7i, O^cj)) -\-Y ( 71 , (j))] 

(C16) 


1 /a 

\ n+1 


^-(n+l) 

2{n + 1)^ Vr 

j [nX{n,e,(p) + Y{n 

■,o,4>)] 




(C17) 


1 / 

n+l 


X-(n+l) 

n{n + 1) \rJ 

1 Z{n,0,,j>), 

(CIS) 


we can obtain the flow field as 


M-Moo = X] ^ (»'X-(n+l)) + V$_(„+i) 

n=l 

n — 2 2 t -7 n + 1 

~7]2n{2n-lf ^“(”+1) + 7?n(2n - 1)’ 

(C19) 
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where 77 is a shear viscosity of the fluid. The viscosity 77 
cancels out and does not appear in the final expression 
of n, while it does appear in the pressure formula jd^ . 

The solved flow field is shown in Fig. [TTJ We plot the 
flow field in the particle frame u and that in the labora¬ 
tory frame u — Uoo- The far-fleld flow in the plane parallel 
to the electric held corresponds to the flow field created 
by pushers Q, because the higher order terms decay 
much faster than the lowest order term representing a 
force dipole. Because the flow in this plane converges in 
y direction, the flow in the plane perpendicular to the 
electric held is divergent. 

We can Anally understand the reason of the instability 
of why the parallel configuration of the two Janus parti¬ 
cles moving in the same direction is unstable (Fig. [12]). 
The hydrodynamic flow created by one particle repels 
the other particle. This interaction makes this configu¬ 
ration unstable and the particles push each other to go 
apart. These kinds of instabilities have been discussed 
especially for the case of pusher-type microswimmers in 
some references [HO, IHli • 

Here we do not take into account the effects of the 
walls or the electrodes, and we cannot explain why the 
particles tend to form clusters moving in the same direc¬ 
tion. Precise calculations or hydrodynamic simulations 
with the walls and the electrostatic interactions might 
account for the cluster formation. 


FIG. 11. (Color online) Streamlines of the flow field around 
a Janus particle. (a)(c)(e) Flow parallel to the electric field. 
(b)(d)(f) Flow perpendicular to the electric field, (a) and 
(b) are the flow field in the particle frame, and (c)-(f) are 
that in the laboratory frame. The far-field flow shown in (e) 
is reminiscent of that of pusher-type microswimmers. The 
radius of the particle is set to 1. The color indicates the 
amplitude of the flow field. The speed of the flow is faster in 
the red region. 

Appendix D: Pair correlation function 

In order to examine the positional order of the crowd of 
the Janus particles, we calculate a pair correlation func¬ 
tion g{r) defined as 



where po is the average number density, N is the total 
number of the particles, and Vij is the distance between 
the i-th and the j-th particles. The brackets mean aver¬ 
age over all the time series. 

The pair correlation obtained from the experimental 
data is similar to that of simple liquids (Fig. [T3|) . It 
means that there is no obvious ordered structure in the 
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FIG. 12. (Color online) The parallel configuration of the 
two Janus particles moving in the same direction is unstable. 
Because the observed plane in the experiment is the plane 
perpendicular to the applied electric held, the how shown in 
this hgure is accountable for the hydrodynamic interaction. 
The hydrodynamic how between the two particles repels each 
other apart. The radius of the particle is set to 1. 


position of the particles like crystals. The prohle of the 


obtained pair correlation is the same as the one observed 
when particles with hnite sizes are randomly distributed. 
This prohle is also observed in the experiments on the 
Quincke rollers in the regime called “polar liquid phase” 
where macroscopically directed motion appeared @ . 
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FIG. 13. (Color online) Pair correlation function g(r) has the 
same form as simple liquids. There is no positional order in 
the crowd of the Janus particles. Error bars are smaller than 
symbols. 
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